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摘要 :【 目的 】 研 究 不 同 茶园 茶 小 绿叶 蝉 Empoasca onukii 种 群 的 遗传 分 化 和 基因 流 格局 ,探究 其 种 
群 遗 传 结构 差异 和 扩散 特点 及 其 与 茶园 景观 格局 的 关系 。【 方 法】 以 福建 省 安溪 县 为 研究 区 域 , 在 
运用 ArcMap10. 5 和 В 软件 包 raster 分 析 研 究 区 域 景观 组 成 和 结构 的 基础 上 ,选取 周边 景观 格局 不 
同 的 18 个 茶园 采集 茶 小 绿叶 蝉 (530 个 个 体 ) ,基于 23 个 微 卫星 位 点 对 这 些 个 体 进行 PCR 扩 增 和 
基因 型 测定 ,采用 ARLEQUIN 3.5.2, FSTAT 2.9.3 fo В 软件 包 adegenet 2.0.0 分 析 其 遗传 多 样 性 
和 遗传 分 化 情况 ;采用 STRUCTURE 2.3.4 fe R 软件 包 adegenet 2.0.0 中 的 DAPC 程序 分 析 其 种 群 
遗传 分 化 ;利用 BAYESASS 3.0.4 估算 种 群 间 最 近 几 代 的 迁移 率 ; 同 时 利用 景观 遗传 学 的 统计 方法 
Mantel, 3E 8j 45 № 3 Am JIRA (MRM) 和 一 般 线性 模型 (GLMM ) 将 茶 小 绿叶 蝉 种 群 遗 传 分 化 与 景 
观 组 成 和 景观 结构 的 空间 数据 进行 关联 分 析 。【 结 果 】 供 试 18 个 茶 小 绿叶 蝉 种 群 的 23 个 微 卫 星 
位 点 的 等 位 基因 数 为 9~52, 等 位 基因 丰富 度 (A) 为 3.686 ~4.397, 基 因 多 样 性 (Di) 为 0.676~ 
0.734, 期 望 杂 合 度 (He) 为 0.659 ~0.729。 聚 类 分 析 显 示 , 有 大 量 森 林 生 境 的 西北 部 5 个 相 邻 样 点 
的 种 群 组 成 一 个 类 群 ,有 大 量 居民 点 和 其 他 作物 田 分 布 的 东北 部 的 5 个 相 邻 样 点 的 种 群 分 为 两 个 
类 群 ,而 集约 化 种 植 模式 明显 的 南部 样 点 的 种 群 也 分 成 两 个 类 群 。 两 两 种 群 间 的 Nei 氏 遗 传 距离 
为 0.042 ~ 0. 984, Provesti K 3t 4 3E 0. 207 ~ 0. 650, 遗传 分 化 指数 sy 为 0. 002 ~ 0.222, 
BAYESASS 分 析 显 示 ,种 群 间 的 现时 基因 流 (Nm) 相对 较 低 , 介 于 0.007 ~0.180, 而 在 种 群 内 部 的 
基因 流 介 于 0. 674 ~0.854。 在 此 基础 上 ,Mantel 和 距离 矩阵 多 元 回归 模型 分 析 显 示 , 地 理 距 离 和 采 
样 点 周围 1 000 m 范 围 内 草地 面积 占 比 是 与 茶 小 绿叶 蝉 种 群 遗传 分 化 有 关 的 两 个 关键 因子 ;一 般 线 
性 模型 分 析 进 一 步 证 实 茶 小 绿叶 蝉 种群 遗 传 分 化 与 这 两 个 关键 因子 线性 关系 显著 , 茶 小 绿叶 蝉 种 
群 遗 传 多 样 性 与 采样 点 周围 半径 1 000 m 范围 内 草地 面积 占 比 线性 关系 显著 。【 结论 】 结 果 说 明 研 
究 区 域内 茶 小 绿叶 蝉 种 群 形成 了 明显 的 遗传 分 化 结构 和 地 理 距 离 隔 离 格局 ,异地 种 群 间 现 时 基因 
交流 受 限 明 显 ,因此 推测 茶 小 绿叶 蝉 其 自主 扩散 能 力 比 较 有 限 , 不 具有 远 距 离 迁 飞 的 习性 。 茶 园 周 
围 1 000 m 范围 内 草地 生境 对 茶 小 绿叶 蝉 种 群 的 遗传 多 样 性 具有 积极 作用 ,茶园 周边 景观 组 成 和 
尺度 范围 可 能 通过 影响 茶 小 绿叶 蝉 种 群 的 扩散 和 定 殖 过 程 而 与 其 种 群 遗 传 结构 相关 联 。 

关键 词 : 茶 小 绿叶 蝉 ; 茶园 ; 景观 格局 ; 遗传 结构 ; 基因 流 ; 害虫 可 持续 治理 
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Crops, Ministry of Agriculture, Fujian Agriculture and Forestry University, Fuzhou 350002, China) 
Abstract: [Aim] This study aims to analyze the patterns of genetic differentiation and gene flow among 
different populations of the tea green leafhopper, Empoasca onukii, in different tea plantations, and to 
investigate the genetic structure difference and dispersion characteristics of E. onukii populations and 
their relationship with the landscape pattern around tea plantation. [Methods] A total of 530 Е. onukii 
samples were collected from 18 tea plantations in Anxi County, Fujian Province, southeastern China. 
ArcMap 10.5 and raster of R package were used to analyze the landscape composition and configuration 
around each studied tea plantation. Е. onukii samples were then individually amplified and individually 
genotyped at 23 microsatellite loci. The genetic diversity and genetic differentiation of E. onukii 
populations were measured using ARLEQUIN 3.5.2, FSTAT 2. 9. 3 and В package adegenet 2. 0. 0. 
Population genetic differentiation was analyzed using both STRUCTURE 2.3.4 and DAPC in R package 
adegenet 2. 0. 0. BAYESASS 3.0.4 was used to estimate the dispersal rates over the last few generations 
among populations. Meanwhile the correlation between the genetic differentiation of E. onukii populations 
and the spatial data of landscape composition and landscape configuration was analyzed by using Mantel 
tests, multiple regression on distance matrices ( MRM) and generalized linear mixed models ( GLMM). 
[Results] All the 23 microsatellite loci in 18 Е. onukii populations were polymorphic, with alleles 
varying from 9 to 52, the allelic richness (Az) from 3. 686 to 4. 397, the gene diversity ( Dıy) from 
0. 676 to 0. 734, and the expected heterozygosity ( He) from 0. 659 to 0. 729. Cluster analysis showed 
that the populations of five sampling sites from the tea plantations surrounded by woodland habitats in the 
northwestern Anxi clustered into one distinct cluster, the populations of five sampling sites from the tea 
plantations surrounded by residential areas and croplands in the northeastern Anxi formed two separate 
clusters, and the populations of sampling sites from the intensified tea plantations in the southern Anxi 
formed two clusters. The Nei’ s distance between two populations ranged from 0. 042 to 0. 984, Provesti' s 
genetic distance from 0. 207 to 0. 650, and the pairwise К.’ value from 0. 002 to 0. 222. BAYESASS 
analysis showed fairly low gene flow (Nm) (0. 007 - 0. 180) among populations, and the gene flow 
within population ranged from 0. 674 to 0. 854. The Mantel tests and MRM model demonstrated that 
geographical distance and the area percentage of grassland in the circular sector of 1 000 m radius around 
sampling site were two key factors in shaping genetic patterns of E. onukii populations. GLMM further 
confirmed a significant linear relationship between the genetic differentiation of E. onukii populations and 
these two key factors, and also between genetic diversity of Е. onukii populations and the area percentage 
of grassland in the circular sector of 1 000 m radius around sampling site. [ Conclusion] The results 
provide direct evidence for significant population genetic differentiation and isolation by distance in E. 
onukii, and suggest that limited recent gene flow has occurred among populations sampled from different 
sites. These findings support the idea that E. onukii may have limited dispersal capacity and is unlikely 
to naturally undergo long-distance migration. The grassland habitat adjacent to tea plantation at the spatial 
scale of 1 000 m has positive effects on the genetic diversity of E. onukii populations. The landscape 
composition and scale around tea plantations may be related to the genetic structure of E. onukii by 
affecting the dispersal and colonization of E. onukii populations. 

Key words: Empoasca onukii; tea plantation; landscape pattern; genetic structure; gene flow; 
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由 于 茶树 为 多 年 生 常 绿 木 本 植物 ,传统 茶园 中 
的 生物 群落 结构 通常 比 其 他 作物 生态 系统 的 群落 结 
构 复 杂 ,并 且 在 相当 长 的 时 间 内 可 以 保持 动态 平衡 
( 陈 李 林 等 , 2019) 。 近 年 来 , 随 着 茶树 种 植 面积 的 
持续 增加 ,茶树 栽培 制度 也 发 生 了 变化 , 越 来 越 多 地 
采用 集约 化 经 营 和 单一 化 种 植 的 模式 ,导致 茶园 及 
其 周围 的 景观 环境 也 逐渐 趋 于 简单 化 ,容易 引发 茶 
园 病虫害 的 流行 和 扩散 ( 李 乐 ,2012; AFT, 
2014; 陈 李 林 等 , 2019 ) 。 与 此 同时 ,大 量 使 用 化 肥 
和 农药 致使 茶园 生态 平衡 日 益 遭 到 破坏 ,引起 有 害 
生物 特别 是 刺 吸 式 害 虫 的 发 生 危 害 不 断 加 

茶 小 绿叶 蝉 Empoasca onukii 是 茶园 中 分 布 最 
广 ,为 害 持续 时 间 最 长 ,造成 损失 也 最 大 的 头号 害 
虫 ,在 我 国 可 造成 减产 10% ~ 1596 ,严重 时 减产 
50% 以 上 (Chen et al., 2019). 。 有 学 者 提出 茶 小 绿 
叶 暗 是 在 20 世纪 60 年 代 侵入 茶园 后 迅速 上 升 为 我 
国 茶 园 害虫 区 系 中 的 头号 害虫 ( ЕЖК, 1979) ,并 
普遍 认为 这 可 能 与 茶园 景观 环境 的 简单 化 过 程 以 及 
化 肥 、 农 药 施 用 水 平 的 变化 有 关 ( 李 乐 , 2012; 周 宁 
T, 2014; Chen et al., 2019)。 因 此 ,在 周边 景观 格 
局 分 布 不 同 的 茶园 中 对 茶 小 绿叶 蝉 种 群 结构 差异 和 
扩散 特点 进行 研究 有 助 于 理解 其 发 生机 制 ,进而 为 
茶园 景观 合理 布局 和 开展 可 持续 害虫 治理 提供 重要 
的 信息 。 

ir 20 年 来 ,得 益 于 分 子 技术 和 地 理 空间 分 析 的 
快速 发 展 , 分 子 水 平 上 的 微观 分 析 手 段 与 景观 生态 
学 的 宏观 统计 工具 相 结 合 ,对 于 连接 景观 环境 与 生 
物种 群 结构 .适应 性 和 迁移 扩散 过 程 提 供 很 大 的 帮 
HJ] ( Storfer et al., 2007) 。 一 些 学 者 利用 RAPD 线 粒 
体 基因 序列 和 微 卫星 分 子 标记 对 茶 小 绿叶 蝉 不 同 地 
理 种 群 进行 了 分 析 ( 李 乐 , 2012; Fu et al., 2014; 周 
TFF, 2014; Zhang et al., 2016, 2019) 。 这 些 研 究 
的 重点 大 多 用 于 解析 不 同 种 群 的 亲缘 关系 和 确认 是 
否 属于 相同 的 优势 种 , 较 少 通 过 基因 流 的 估算 和 基 
因 流 格局 的 分 析 以 获得 其 种 群 迁 移 、 扩 散 相 关 的 信 
息 , 也 没有 将 这 些 遗 传 结构 基因 流 和 扩散 信息 与 茶 
园 景观 环境 相关 联 , 并 用 于 阐明 茶 小 绿叶 蝉 在 不 同 
茶园 景观 环境 中 的 种 群 分 布 和 扩散 特点 。 本 研究 以 
福建 省 安溪 县 为 研究 区 域 ,选取 了 周边 景观 结构 和 
组 成 不 同 的 18 个 茶园 采集 茶 小 绿叶 蝉 ,基于 23 个 
微 卫 星 分 子 标记 对 这 些 茶 小 绿叶 蝉 种 群 进行 了 遗传 
结构 和 基因 流 格局 的 分 析 , 探 究 在 周边 景观 格局 不 
同 的 茶园 中 茶 小 绿叶 暗 的 种 群 结构 差异 和 扩散 特 
点 ,为 阐明 其 灾变 机 理 \ 优 化 茶园 景观 合理 布局 方案 
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和 害虫 可 持续 治理 提供 依据 。 
1 材料 与 方法 


1.1 研究 区 域 选择 和 样本 采集 

选择 福建 省 安溪 县 作为 研究 区 域 ,县域 范围 
117°36' ~ 118°17'Е, 24°50’ ~ 25"26'N ,全 县 总 面积 
3 057.28 km ( 税 伟 等 , 2017) 。 区 域内 地 势 从 西北 
向 东南 倾斜 ,西北 部 山峰 林立 ,东南 部 地 势 相 对 平 
缓 ,以 丘陵 山地 为 主 。 茶 园 在 全 县 区 域内 均 有 分 布 ， 
但 由 于 明显 的 地 势 地 貌 差 异 ,不 同乡 镇 的 种 植 模 式 
有 明显 差异 :在 位 于 南部 的 西 坪 等 乡镇 ,主要 采用 单 
一 化 和 集约 化 的 种 植 模式 ,在 茶园 周围 少见 自然 和 
半 自 然 生境 或 其 他 作物 田 ;在 位 于 西北 部 的 桃 入 和 
福田 等 乡镇 ,茶园 大 都 分 布 于 高 海拔 的 山坡 上 ,周围 
经 常会 有 和 森林 等 自然 和 半 自 然 生境 分 布 ;而 在 位 于 
东部 的 蓬莱 和 湖 头 等 乡镇 ,种 植 模式 则 会 更 多 样 化 ， 
茶园 周围 经 常会 有 其 他 多 种 作物 田 的 分 布 ( 税 伟 
等 , 2017) (#1). 

基于 以 上 信息 ,在 研究 区 域内 选取 了 周边 景观 
格局 不 同 的 18 个 茶园 (其 中 茶树 有 多 年 树龄 ) 进行 
研究 ( 表 1) , 避 开 选择 新 建 茶园 以 避免 近期 茶 苗 移 
植 对 分 析 的 影响 。2018 4E 8 月 ,在 选 定 的 每 个 茶园 
中 ,随机 选择 位 于 其 中 央 的 5 ~ 10 个 点 (彼此 间隔 
20 ~ 100 m) ,用 扫 网 法 采集 至 少 50 KAR SRM Л, 
虫 。 本 研究 将 来 自 同一 样 点 的 个 体 作为 一 个 种 群 用 
于 后 续 的 研究 和 分 析 。 

1.2 微 卫 星 位 点 基因 型 检测 

本 研究 使 用 天 根 生化 科技 (北京 ) 有 限 公 司 的 
微量 样品 基因 组 提取 试剂 盒 , 对 单 头 茶 小 绿叶 蝉 进 
行 基因 组 DNA 的 提取 。 为 了 保证 样品 种 类 的 一 致 
性 ,分析 个 体 在 基因 组 提取 之 前 首先 进行 初步 形态 
学 鉴定 (Qin et al., 2015), 并 参照 Demichelis 等 
(2010) 的 СОТ 基因 引物 以 及 PCR 反应 体系 和 条 件 
对 每 份 基因 组 样品 进行 扩 增 ,将 产物 序列 (经 ABI 
3730XL 测序 成 功 ) Ej NCBI 中 茶 小 绿叶 蝉 Empoasca 
onukii 的 线粒体 全 长 序列 (Liu et al., 2017) 进行 比 
对 , 仅 保留 可 以 比 对 上 的 样品 用 于 后 续 微 卫星 位 点 
的 扩 增 和 分 析 。 

本 研究 选取 23 个 微 卫星 位 点 作为 分 子 标记 研 
究 茶 小 绿叶 蝉 种 群 遗 传 变异 ,引物 序列 信息 参照 
Zhang 等 (2016) 和 李 傅 等 (2016)( 表 2)。PCR 扩 增 
采用 天 根 生化 科技 (北京 ) 有 限 公 司 的 多 重 PCR 扩 
增 试剂 盒 (KT109) ,将 4 ~6 对 引物 组 合 在 同一 个 体 
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A1 福建 省 安溪 县 茶 小 绿叶 蝉 采样 点 信息 
Table 1 Information of sampling sites of Empoasca onukii in Anxi County, Fujian Province, southeasterm China 
种 群 样本 采集 地 经 纬度 海拔 (m) 采集 日 期 (年 .月 ) 
Population Sampling site Longitude and latitude Altitude Collection date (year. month ) 
FX 湖上 Hushang 25.25°N, 118.00°E 530 2018.8 
MB ЗЕЕ Penglai 25. 14?N, 118. 11°Е 136 2018.8 
SL 金谷 Jingu 25.21°N, 118. 10E 189 2018.8 
SG Ей Bailai 25.26°N, 118.05°E 439 2018.8 
HL М Shangqing 25.18°N, 117.98°Е 785 2018.8 
XH 祥 华 Xianghua 25.18?N, 117. 77?E 857 2018.8 
HM 长 坑 Changkeng 25.18°N, 117. 84°E 695 2018.8 
WZ 监 田 Lantian 25.09°N, 117. 88°Е 502 2018.8 
HZ 感 德 Gande 25.27°N, 117.85°Е 542 2018.8 
XG BERE Taozhou 25.31°N, 117.81?E 925 2018.8 
JD 剑 斗 Jiandou 25.32°N, 117.95°Е 300 2018.8 
FQ 福田 Futian 25.26°N, 117. 65°Е 540 2018.8 
JB 6 Ниси 25.02?N, 118. 00°Е 286 2018.8 
ҮҮ QE Xiping 25.01?N, 117. 94°Е 607 2018.8 
ІТ P5 FA Lutian 25.03?N, 117. 89°Е 502 2018.8 
MS KLE Daping 24.93°N, 117.96°Е 963 2018.8 
XD 龙门 Longmen 24.93°N, 118. 02°Е 710 2018.8 
JY JEiH Longjuan 24.96°N, 117. 790E 565 2018.8 
表 2 本 研究 中 所 用 的 23 个 微 卫星 位 点 的 引物 信息 
Table 2 Primer information of 23 microsatellite loci in this study 
l " Е 荧光 标记 "m 
位 点 重复 碱 基 引物 序列 〈5 -37) 产物 大 小 (pp) 参考 文献 
Florescent 
Loci Repeat motif Primer sequence Product size References 
marker 
F: TACTCACCGCTCCTCTATC 
Eo-1-57 (АС), 93 – 105 5'-ТАМКА Zhang et al., 2016 
В: CACTTTTATTTTCGGCTCT 
F: GTTTTGGCTATCATTAG 
Eo-37 (АС), 159 – 169 5'-НЕХ Zhang et al., 2016 
В: ССАСТСАССАААССАС 
Е; ССААССССТТАССАТАСТ 
Eo-E-12 (6T), 123 -165 5'-TAMRA Zhang et al., 2016 
R: AACACCGCATTTCATACA 
F: СССААСААСТССААСАТАА 
Eo-1-5 (GTT), 121 -145 5'6-FAM Zhang et al., 2016 
R: CGAGGAGGAGGCTAAAG 
F: CGGCATTCATTATCTC 
Eo-1-61 (CA); (CA), 84 -106 5'-HEX Zhang et al., 2016 
R: CACACAACTCACTCGCT 
F: GCCGTGTGTAATGGTATCC 
Eo-1-52 (AC) i 197 -217 5'-HEX Zhang et al., 2016 
R: GACGCCTAGCAATGTT 
F: GGTTTGAATGGGCTCAGATT 
EVMS10 (AAT); 145 5'-TAMRA 李 倩 等 , 2016 
В: CCTAGGAAAGCGTGGTCAAG 
Е; AACACATCACATGCCACCAT 
EVMSOI (CT); 2778 5'-TAMRA 李 倩 等 , 2016 
В: CCTAAATGTAACCATTCACCTGC 
F: GATGGAGAAGGAATGTTCGC 
EVMS02 (GCT); 213 5'6-FAM 李 倩 等 , 2016 
В: СТССАСАААССАТСААСССА 
Е; ССАСАСАСССААСТСААТ 
Eo-36 (AG), É 265 -295 5'-ROX Zhang et al., 2016 


: GGTACCCGAAGGAAGGAT 
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续 表 2 Table 2 continued 
Я 2: zs 荧光 标记 
位 点 重复 碱 基 引物 序列 (5' -3 ) 产物 大 小 (bp) 参考 文献 
Florescent 
Loci Repeat motif Primer sequence Product size References 
marker 

F: CAGGAGCAGGACAAGAG 

Eo-70 (AC), C(CA) 12 211 -257 5'-TAMRA Zhang et al., 2016 
R: GCACATAAGCCTAAACAGAC 
F: GATTTTTGGGCTCCTGTCAA | 

ЕУМ505 (TAA); 158 5'-HEX 李 倩 等 , 2016 
R: CATGTATTATTAGACTTGGTCAAGCA 
Е; CTTCGAGTCCGCTTTGAAAC 

EVMS06 (TGA); 108 5'6-FAM 李 倩 等 , 2016 
В: TGCTCATGCGTAATTGCTTC 
F: ССАСТСАССАААССАССА 

Ео-29 (TC), 154 - 168 5'6-FAM Zhang et al., 2016 
R: GGGTATGATTAGCGGTGT 
F: TCGCTCCACTCTACCACT 

Eo-51 (CT), GT(CT), 194 -220 5'-ROX Zhang et al., 2016 
В: СТААААСТААСАААТССАССТ 
F: СССАТААССААССТААСС 

Ео-9 (ТС), 96 – 128 5'6-FAM Zhang et al., 2016 
R: AGGAGTCGGGAGTCTGTT 
Е; TAACAGATGAGGGTAGGAAT 

Eo-1-65 (CA), 208 -228 5'-HEX Zhang et al., 2016 
В: CTTTAGAAATGGCAAGACAG 
F: TGACAGCCATAAACACCG 

Eo-1-82 (GT); 274 – 304 5'-ROX Zhang et al., 2016 
В: CGTAGACCAGATGACCCTC 
F: CTGTCACTTGCCAATAACTCT 

Eo-20 (GA), 125 - 155 5'-ROX Zhang et al., 2016 
R: CAACCACCTCACTCCCTCT 
F: CTGTTCGCAGTTCACATCATTC 

Eo-54 (СА), 403 - 443 5'6-FAM Zhang et al., 2016 
R: GACCCGCTACGCTTACCTATT 
F: TAGGGGTTCGACAGACTTG 

Eo-68 (AG) ,GG(AG), 154 - 192 5'-TAMRA Zhang et al., 2016 
R: GAGGTGAGGTTGGGATTTG 
F: GAGTCTGTTCGGTTTGATGT 

Eo-83 (АС), 84 – 108 5'-ROX Zhang et al., 2016 
В: СААСССАТААССААССТААС 
F: GGACAATGAAAATTCGAGGACGG 

Eo42 (GA) g 117 - 141 5'-HEX Zhang et al., 2016 
R: GGCAATCGGCAACAACAAAC 


系 中 进行 扩 增 ,反应 体系 和 条 件 参 照 其 使 用 说 明 。 
PCR 产物 经 电泳 检测 后 ,利用 ABI 3730XL 测序 仪 
( Applied Biosystems ) 结合 GENEMAPPER 3. 2 软件 
( Applied Biosystems) 分 析 确 定 其 中 对 应 位 点 的 片段 
长 度 和 识别 片段 基因 型 ,利用 Micro-Checker ( Van 
Oosterhout et al., 2004) 检测 是 否 存在 因为 无 效 等 位 
基因 和 滑 移 导 致 的 位 点 读 取 错 误 的 问题 。 最 终 ,我 
们 获得 了 来 自 安溪 县 18 个 样 点 的 530 个 个 体 的 微 
卫星 基因 型 数据 用 于 后 续 种 群 分 析 。 
1.3 景观 格局 分 析 

为 了 定量 化 解析 茶园 周边 景观 格局 特征 与 茶 小 
绿叶 蝉 种 群 遗 传 信息 之 间 的 相关 性 ,从 地 理 数据 云 
EG (Һир: /www. gscloud. cn/ ) 下 载 获 取 了 2018 年 
国家 标准 土地 利用 类 型 分 类 (land cover) 原始 影像 
图 ,其 分 辨 率 大 约 为 30 m, 包 括 了 6 大 类 国家 标准 
的 土地 利用 类 型 分 类 。 在 该 影像 图 的 基础 上 ,结合 
谷歌 地 球 ( Google Earth ) 的 投影 影像 和 采样 实地 考 


























察 的 照片 资料 ,把 研究 区 域内 茶园 单独 色 绘 出 来 , 构 
建 了 包含 茶园 .林地 、 草 地、 耕地 水域 .建筑 和 未 利 
用 ( 裸 地 、 岩 石 等 )7 个 土地 利用 类 型 的 矢量 影像 资 
料 。 利 用 该 矢量 数据 ,我 们 对 研究 区 域 的 景观 组 成 
和 结构 进行 了 分 析 。 

使 用 ArcMapl0.5 分 别 生 成 统计 了 各 采样 点 周 
围 半 径 500, 1 000, 1 500, 2 000 和 2 500 m 内 各 类 
土地 利用 类 型 的 面积 占 比 (分 析 中 将 水 域 . 建 筑 和 
未 利用 合并 为 一 个 无 植被 的 土地 利用 类 型 统计 )， 
使 用 R(R Core Team, 2018) 中 的 cor 函数 对 25 个 
表示 景观 组 成 的 变量 进行 了 相关 性 检验 ,只 保留 没 
有 明显 共 线性 关系 (r<0.8) 的 11 个 景观 组 成 因子 
( 表 3)。 同 时 ,使 用 R 软件 包 raster(Hijmans and 
van Etten, 2016) 将 景观 空间 数据 转化 生成 阻力 距 
离 曲 面 的 栅 格 数据 ,以 表示 不 同 土地 利用 类 型 对 茶 
小 绿叶 蝉 扩散 的 阻力 大 小 (Chen et al., 2019) ,分 析 
获得 代表 景观 结构 的 因子 ( 表 3 ) 。 此 外 ,使 用 手持 
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RI 本 研究 中 所 用 的 景观 环境 变量 信息 
Table 3 Landscape and environmental variables used in this study 
变量 变量 描述 类 别 及 意义 
Variable Description of variable Category 
地 理 距离 
geo 地 理 距离 Geographical distance 
Geographical distance 
А 采样 点 周围 半径 500 m 内 茶园 面积 占 比 景观 组 成 
500m_Tea plantation 
Area percentage of tea plantations in the circular sector of 500 m radius around sampling site Landscape composition 
采样 点 周围 半径 500 m 内 耕地 面积 占 比 景观 组 成 
500m_Cropland Я 
Area percentage of cropland in the circular sector of 500 т radius around sampling site Landscape composition 
采样 点 周围 半径 500 m 内 草地 面积 占 比 景观 组 成 
500m_Grassland 
Area percentage of grassland in the circular sector of 500 m radius around sampling site Landscape composition 
采样 点 周围 半径 500 m 内 森林 面积 占 比 景观 组 成 
500m_Woodland 
Area percentage of woodland in the circular sector of 500 m radius around sampling site Landscape composition 
采样 点 周围 半径 1 000 m 内 草地 面积 占 比 景观 组 成 
1000m_Grassland 
Area percentage of grassland in the circular sector of 1 000 m radius around sampling site Landscape composition 
采样 点 周围 半径 1 500 m 内 森林 面积 占 比 景观 组 成 
1500m_Woodland 
Area percentage of woodland in the circular sector of 1 500 m radius around sampling site Landscape composition 
采样 点 周围 半径 2 500 m 内 茶园 面积 占 比 景观 组 成 
2500m_Tea plantation 
Area percentage of tea plantations in the circular sector of 2 500 m radius around sampling site Landscape composition 
采样 点 周围 半径 500 m 内 无 植被 覆盖 土地 面积 占 比 
景观 组 成 
2500m_noveg Area percentage of artificial and open spaces without vegetation in the circular sector of 2 500 
Landscape composition 
m radius around sampling site 
采样 点 周围 半径 2 500 m 内 耕地 面积 占 比 景观 组 成 
2500m_Cropland 
Area percentage of cropland in the circular sector of 2 500 m radius around sampling site Landscape composition 
采样 点 周围 半径 2 500 m 内 草地 面积 占 比 景观 组 成 
2500m_Grassland 
Area percentage of grassland in the circular sector of 2 500 m radius around sampling site Landscape composition 
采样 点 周围 半径 2 500 m 内 森林 面积 占 比 景观 组 成 
2500m_Woodland 
Area percentage of woodland in the circular sector of 2 500 m radius around sampling site Landscape composition 
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Landscape configuration 


景观 结构 


Landscape configuration 


景观 结构 


Landscape configuration 


В 软件 包 geosphere 1. 5-10 ( Hijmans et al., 2017) it 


高 度 和 经 纬度 坐标 ( 表 1) ,并 基于 经 纬度 信息 使 用 算 采 样 点 间 地 理 距 离 , 从 而 获得 海拔 和 地 理 距离 的 
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因子 ( 表 3) 。 
1.4 种 群 遗传 多 样 性 分 析 

利用 ARLEQUIN 3. 5. 2( Excoffier and Lischer, 
2010) ib Я © Th HE HH E Ze d HE ( expected 
heterozygosity, Не) #1 1 2 Z& £r JE ( observed 
heterozygosity, Ho) ,用 FSTAT 2.9.3( Goudet, 2001) 
计算 各 种 群 的 基因 多 样 性 ( gene diversity, Dy) 和 等 
位 基因 丰富 度 (allelic richness, Ах) 等 遗传 多 样 性 
参数 。 
1.5 ”种群 遗传 结构 分 析 

使 用 STRUCTURE 2. 3. 4 ( Pritchard et. al., 
2000) xf 18 个 种 群 进行 贝 叶 斯 聚 类 分 析 。 选 择 R 
软件 包 adegenet 2. 0.0(Jombart，2008 ) 中 不 基于 模 
型 假设 的 DAPC ( discriminant analysis of principal 
components ) (Jombart et al., 2010) 程序 解析 个 体 之 
间 的 遗传 关系 。 
1.6 ”种群 遗传 分 化 分 析 

利用 В 软件 包 adegenet 2.0.0 中 的 dist. genpop 
函数 计算 种 群 间 的 遗传 距离 ,包括 Nei 氏 遗 传 距离 
(Nei' s genetic distance) ( Nei, 1972, 1978) 和 
Provesti 氏 遗 传 距离 ( Provesti' s genetic distance ) ,并 
用 ARLEQUIN 3. 5. 2( Excoffier and Lischer, 2010) if 
算 种 群 的 遗传 分 化 指数 (pairwise For) 和 显著 性 检 












































(multiple regression on distance matrices, MRM) 分 析 
(Wang, 2013) ,以 进一步 确认 与 研究 区 域内 茶 小 绿 
叶 蝉 种群 Provesti 氏 遗 传 距 离 、Nei 氏 遗 传 距离 和 遗 
传 分 化 指数 有 关 的 景观 环境 因子 。 

使 用 R 软件 包 stats 的 一 般 线 性 模型 
( generalized linear mixed model, GLMM) 分 析 采 样 点 
周围 半径 1 000 m 范围 内 的 草地 面积 占 比 (1000m_ 
Grassland) 和 地 理 距离 (geo) 与共 小 绿叶 蝉 种 群 
Provesti 氏 遗 传 距离 、Nei 氏 遗 传 距离 和 遗传 分 化 指 
数 的 线性 关系 ,并 分 析 各 样 点 1000m_Grassland 值 与 
茶 小 绿叶 蝉 种 群 遗传 多 样 性 ( Dj,, Ак, He, Ho) 的 线 
性 关系 。 

















2 结果 


2.1 种 群 遗传 多 样 性 

23 个 微 卫星 位 点 在 18 个 种 群 中 均 表 现 出 多 态 
性 ,等 位 基因 数量 从 9(Eo-1-57) 到 52( Eo-51) 8; 
等 位 基因 丰富 度 (hx) 介 于 3. 686 (FQ) 和 4. 397 
(HL) 之 间 , 平 均值 为 4.5195 ;基因 多 样 性 (Di) 介 
于 0.676(FX) 和 0.734(HL) ,平均 值 为 0.7081; 期 
BAAR GE (He) Jt 0. 659 (FQ) F 0. 729 ( HL) ,平均 
值 为 0.7567( 表 4) ,位 于 安溪 县 西部 和 北部 样 点 的 








Jy PH. 
1.7 种 群 基因 流 格局 分 析 

使 用 BAYESASS 3. 0. 4 ( Wilson and Rannala, 
2003 ) 估算 样 点 之 间 茶 小 绿叶 蝉 种 群 最 近 几 代 的 迁 
移 率 (m) , 即 每 个 世代 中 迁移 个 体 在 种 群 中 所 占 的 
比例 。 设 置 程序 运行 10 次 ,每 次 进行 20 000 000 次 
重复 的 MCMC 搜索 ( MCMC interaction) ,舍弃 最 初 
的 2 000 000 次 重复 的 数据 ,每 隔 100 次 采样 一 次 以 








确保 m 值 的 可 靠 \ 准 确 。 
1.8 茶园 周边 景观 格局 与 茶 小 绿叶 蝉 种 群 遗传 的 
相关 性 分 析 


以 11 个 景观 组 成 因子 和 3 个 景观 结构 因子 以 
及 海拔 (elevation ) 和 地 理 距 离 (geo)( 表 3) 为 解释 
变量 ( explanatory variable) ,使 用 В 软件 包 ecodist 
(Goslee and Urban, 2007) 分 析 其 与 茶 小 绿叶 蝉 种 群 
的 Provesti 氏 遗 传 距离 Nei 氏 遗 传 距离 和 和 遗传 分 化 
指数 之 间 的 Mantel 相关 性 。 

VA 16 个 景观 环境 因子 ( 表 3 ) 为 解释 变量 , 以 
Provesti 氏 遗 传 距 离 Nei 氏 遗 传 距 离 和 遗传 分 化 值 
为 响应 变量 (response variable) ,使 用 R 软件 包 
ecodist 的 MRM 函数 进行 距离 矩阵 多 元 回归 











茶 小 绿叶 蝉 种 群 (如 FQ, XG, HZ 等 ) 表 现 出 较 低 的 
遗传 多 样 性 水 平 。 
2.2 ”种群 遗传 结构 

STRUCTURE 软件 推断 的 天 的 最 大 可 能 值 是 2， 
表明 研究 区 域内 茶 小 绿叶 蝉 种 群 可 以 划分 为 两 个 遗 
传 类 群 (cluster) :一 个 类 群 包括 安溪 县 北部 10 个 样 
点 的 种 群 , 男 一 个 类 群 包括 安溪 县 南部 8 个 样 点 的 
种 群 (图 1: A, F)。 在 此 基础 上 ,K=3 ~6 时 ,安溪 
县 茶 小 绿叶 蝉 种 群 遗 传 分 化 格局 比较 明显 (图 1)， 
大 体 上 为 位 于 安溪 县 西北 部 5 个 相 邻 样 点 (FQ， 
XG, HZ, XH 和 HM) 的 茶 小 绿叶 暗 种 群 组 成 一 个 
类 群 , 分 别 位 于 北部 和 中 部 的 JD 和 WZ 种 群 各 自 单 
独 成 一 个 类 群 ,位 于 东部 的 5 个 相 邻 样 点 的 种 群 
(FX, SG, HL, SL 和 MB) 组 成 一 个 类 群 ,位 于 南部 
的 3 个 相 邻 样 点 的 种 群 (LT, YY 和 MS) 组 成 一 个 
类 群 ,同样 位 于 南部 的 另外 3 个 不 相 邻 样 点 的 种 群 
(JY, 卫 和 XD) 组 成 一 个 类 群 。 结 合 安溪 县 土地 利 
用 分 类 图 来 看 ,位 于 东部 的 5 个 样 点 (FX, SG, HL, 
SL fll MB) 以 及 FQ 和 WZ 周围 的 茶叶 种 植 较 少 ,有 
较 多 的 其 他 土地 利用 类 型 的 分 布 , 与 其 他 样 点 之 间 
处 于 相对 隔离 的 状态 ( 即 与 其 他 样 点 间 没 有 适宜 生 
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表 4 


福建 安溪 茶 小 绿叶 蝉 种 群 的 遗传 多 样 性 


Table 4 Genetic diversity of Empoasca onukii populations in Anxi, Fujian 











种 群 个 体 数 平均 等 位 基因 数 基因 多 样 性 等 位 基因 丰富 度 期 望 杂 合 度 
Population Number of individuals Average number of alleles Gene diversity Dy Allelic richness Аъ Expected heterozygosity He 
FX 32 10.217 0.676 4.054 0.672 
MB 24 9.522 0.716 4.269 0.710 
SL 32 9.087 0. 682 4.010 0. 678 
SG 32 10.174 0.712 4.187 0.708 
HL 32 11.000 0.734 4.397 0.729 
XH 32 10.522 0.705 4.272 0.701 
HM 32 10.435 0.714 4.256 0. 709 
WZ 21 8.739 0.711 4.116 0.704 
HZ 32 10.783 0.723 4.240 0.718 
XG 32 10. 565 0.702 4.230 0.697 
JD 32 10.435 0.698 4.178 0.694 
FQ 5 4.000 0.693 3.686 0.659 

JB 32 10.043 0.709 4.137 0. 704 
YY 32 11.478 0.733 4.374 0.728 
LT 32 10. 783 0.720 4.305 0.716 
MS 32 10.522 0.709 4. 190 0. 705 
XD 32 9.913 0.703 4.077 0.698 
JY 32 9.304 0.706 4. 021 0. 701 

















种 群 信 息 见 表 1, 下 同 。For population information, see Table 1. The same below. 





境 茶园 的 连接 )。 据 此 ,我 们 推测 周边 景观 格局 不 
同 的 茶园 中 的 茶 小 绿叶 蝉 种 群 遗传 结构 和 基因 交流 
模式 可 能 存在 差异 。 

不 基于 模型 假设 的 DAPC 分 析 将 安溪 县 18 个 
茶 小 绿叶 蝉 种 群 划 分 为 4 个 类 群 (图 2) :JY, XD 和 
JB 为 一 个 类 群 ,FX,， XG 和 XH 为 一 个 类 群 , FQ, 
MB, WZ, JD 和 HZ 为 一 个 类 群 ,YY, MS, LT, SG 
和 SL 为 一 个 类 群 。 总 体 上 仍然 是 南北 种 群 分 化 
明显 
2.3 












































种 群 遗传 分 化 
两 个 不 同 种 群 间 的 № 氏 遗 传 距离 为 0. 042 
(MS-LT) ~ 0. 984 ( SL-FQ) , Provesti 氏 遗 传 距离 为 
0.207( MS-LT) ~ 0. 650( SL-FQ) (#25). Week 
虽 数 为 0. 002 ( MS-LT) ~ 0. 222(SL-FQ) ,大 部 分 均 
达到 了 统计 学 上 的 显著 性 水 平 (P<0.05)( 表 6)。 
根据 Wright 等 (1978 ) 提出 的 标准 (0 < Fs «0.05, 
种 群 间 分 化 很 弱 ;0. 05 < For «0. 15 ,中 等 分 化 ;0. 15 < 
Fs, <0. 25 ,分 化 大 ;sw >0.25 ,分 化 极 大 ) ,基于 种 
群 遗传 距离 和 遗传 分 化 指数 ,大 部 分 茶 小 绿叶 蝉 种 
群 的 遗传 分 化 表现 出 中 等 到 较 高 的 水 平 。FQ 与 其 
他 种 群 的 群 遗 传 分 化 较 大 ,HM, SL, SG, MS, LT 和 
YY 种 群 间 遗 传 分 化 也 比较 大 。 



































2.4 种 群 基因 流 格 局 

BAYESASS 的 多 次 单独 运行 获得 了 较为 一 致 的 
结果 ,在 检测 的 306 对 异地 种 群 间 ,大 部 分 种 群 间 的 
现时 基因 流 (Nm) 很 小 , 介 于 0.007 ~0.010 之 间 
(228 对 ) ,其 余 70 对 介 于 0.010 ~0.097 之 间 ,8 对 
介 于 0.108 ~0.180 之 间 ; 而 在 种 群 内 部 个 体 之 间 基 
因 流 基本 都 很 大 , 介 于 0.674 ~0.854 之 间 ( 表 7)。 
说 明 茶 小 绿叶 蝉 异地 种 群 之 间 的 现时 基因 交流 可 能 
受 限 ,其 现时 基因 流 主要 发 生 在 种 群 内 部 个 体 之 间 。 
2.5 茶园 周边 景观 格局 与 茶 小 绿叶 蝉 种 群 遗 传 变 
异 的 相关 性 

Simple Mantel 分 析 显 示 地 理 距 离 ( geo) 与 
Provesti 氏 遗 传 距离 (r 20. 2938, P =0.0088)、Nei 
氏 遗 传 距离 (r 20.2015, P=0.0455) 以 及 遗传 分 化 

EA (r20.2052, P=0.0301) 之 间 均 存在 显著 的 相 

关 性 ( 表 8)。 在 排除 了 地 理 距 离 (geo) 的 影响 后 ， 
Partial Mantel 分 析 显 示 ,11 个 景观 组 成 因子 均 与 
Provesti 氏 遗 传 距离 Nei 氏 遗 传 距离 和 遗传 分 化 指 
数 显著 相关 , 而 海拔 (elevation) 和 与 景观 结构 有 关 
的 所 有 因子 均 没 有 检测 到 与 Provesti 氏 遗 传 距离 、 
Nei 氏 遗 传 距离 和 遗传 分 化 指数 有 显著 相关 性 。 
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主 成 分 分 析 特 征 值 
PCA eigenvalues 


判别 分 析 特 征 值 


DA eigenvalues 
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线性 判别 2 Linear discriminant 2 (LD2) 











线性 判别 1 Linear discriminant 1 (LD1) 


2 DAPC 分 析 茶 小 绿叶 蝉 种 群 的 遗传 关系 
Fig. 2 Genetic correlation analysis of Empoasca onukii populations using DAPC 
点 代表 每 个 分 析 个 体 , 椭 圆 代表 种 群 ,椭圆 中 心 点 之 间距 离 的 大 小 代表 种 群 亲缘 关系 的 远近 。Dots represent the study individuals, inertia ellipses 


represent the populations, and the distance between the center points of the ellipses represents the genetic distance between populations. 





表 5 茶 小 绿叶 蝉 不 同 种 群 之 间 的 Nei 氏 遗 传 距离 (下 三 角 ) 和 Provesti 氏 遗 传 距离 (上 三 角 ) 值 
Table 5 Nei’s genetic distance ( below diagonal) and Provesti' s genetic distance ( above diagonal) 
values among different populations of Empoasca onukii 





И ЕХ MB SL SG HL XH HM WZ HZ XG JD FQ JB YY LT MS XD JY 
FX 0.394 0.371 0.378 0.356 0.246 0.363 0.397 0.337 0.236 0.322 0.532 0.350 0.348 0.367 0.339 0.355 0.340 
MB 0.300 0.526 0.491 0.258 0.387 0.459 0.296 0.358 0.367 0.368 0.448 0.393 0.422 0.449 0.436 0.398 0.405 
SL 0.259 0.664 0.275 0.496 0.395 0.286 0.535 0.422 0.403 0.425 0.650 0.402 0.305 0.312 0.316 0.398 0.405 
SG 0.298 0.614 0.091 0.447 0.392 0.243 0.509 0.407 0.386 0.421 0.640 0.372 0.286 0.310 0.302 0.435 0.421 
HL 0.250 0.066 0.555 0.507 0.353 0.436 0.284 0.336 0.342 0.320 0.467 0.374 0.380 0.413 0.410 0.388 0.376 
XH 0.059 0.323 0.279 0.319 0.245 0.367 0.373 0.300 0.265 0.327 0.533 0.374 0.349 0.369 0.345 0.340 0.357 
HM 0.248 0.510 0.096 0.060 0.434 0.257 0.478 0.397 0.384 0.413 0.619 0.358 0.274 0.287 0.264 0.377 0.375 
WZ 0.314 0.093 0.645 0.606 0.115 0.306 0.533 0.348 0.377 0.404 0.425 0.417 0.444 0.461 0.437 0.399 0.398 
HZ 0.141 0.245 0.347 0.344 0.198 0.125 0.305 0.208 0.352 0.294 0.529 0.411 0.348 0.410 0.385 0.391 0.406 
XG 0.045 0.282 0.285 0.307 0.238 0.072 0.265 0.304 0.169 0.354 0.549 0.362 0.342 0.371 0.355 0.351 0.364 
JD 0.146 0.270 0.335 0.363 0.210 0.149 0.308 0.303 0.101 0.192 0.552 0.433 0.381 0.411 0.402 0.433 0.453 
FQ 0.536 0.236 0.984 0.927 0.279 0.548 0.848 0.201 0.436 0.541 0.546 0.549 0.593 0.592 0.593 0.537 0.536 
JB 0.202 0.265 0.282 0.253 0.223 0.227 0.221 0.295 0.281 0.187 0.378 0.484 0.340 0.351 0.362 0.261 0.241 
YY 0.199 0.411 0.100 0.106 0.334 0.214 0.092 0.433 0.226 0.206 0.263 0.708 0.177 0.248 0.255 0.359 0.380 
LT 0.230 0.445 0.114 0.141 0.367 0.249 0.100 0.459 0.304 0.243 0.316 0.727 0.217 0.070 0.207 0.387 0.378 
MS 0.208 0.442 0.115 0.128 0.371 0.221 0.092 0.441 0.276 0.222 0.303 0.724 0.220 0.069 0.042 0.363 0.377 
XD 0.206 0.279 0.372 0.318 0.222 0.191 0.278 0.297 0.271 0.179 0.381 0.468 0.072 0.239 0.290 0.278 0. 248 
JY 0.184 0.270 0.351 0.295 0.225 0.211 0.288 0.263 0.256 0.180 0.380 0.434 0.069 0.239 0.259 0.262 0.075 
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表 6 茶 小 绿叶 蝉 不 同 种 群 之 间 的 遗传 分 化 指数 值 ( 下 三 角 ) 及 其 统计 检验 值 ( 上 三 角 ) 
Table 6 Pairwaise Fs values ( below diagonal) and their statistic P values ( above diagonal) 
among different populations of Empoasca onukii 


种 群 


Populations 


FX MB SL SG HL XH HM WZ HZ XG JD FQ JB YY LT MS XD JY 





FX 0.000 0.000 0.000 0.000 0.126 0.000 0.000 0.000 0.613 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
MB 0.093 0.000 0.000 0.468 0.000 0.000 0.072 0.000 0.000 0.000 0.027 0.000 0.000 0.000 0.000 0.000 0.000 
SL 0.089 0.171 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
SG 0.095 0.151 0.024 0.000 0.000 0.297 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
HL 0.077 0.008 0.147 0.128 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
XH 0.011 0.092 0.090 0.093 0.069 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
HM 0.080 0.132 0.026 0.009 0.113 0.076 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
WZ 0.096 0.016 0.168 0.149 0.025 0.086 0.136 0.000 0.000 0.000 0.243 0.000 0.000 0.000 0.000 0.000 0.000 
HZ 0.043 0.068 0.105 0.096 0.053 0.033 0.085 0.055 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
XG 0.006 0.083 0.092 0.091 0.069 0.015 0.080 0.087 0.049 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
JD 0.047 0.080 0.106 0.106 0.061 0.044 0.091 0.088 0.026 0.060 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
FQ 0.149 0.046 0.222 0.201 0.060 0.139 0.191 0.036 0.107 0.140 0.141 0.000 0.000 0.000 0.000 0.000 0.000 
JB 0.066 0.076 0.090 0.075 0.063 0.068 0.066 0.083 0.081 0.057 0.111 0.123 0.000 0.000 0.000 0.027 0.072 
YY 0.062 0.107 0.027 0.025 0.087 0.061 0.019 0.110 0.061 0.059 0.076 0.160 0.049 0.054 0.036 0.000 0.000 
LT 0.073 0.117 0.033 0.039 0.097 0.072 0.024 0.119 0.084 0.072 0.092 0.168 0.063 0.011 0.874 0.000 0.000 
MS 0.069 0.120 0.034 0.035 0.101 0.067 0.021 0.119 0.079 0.068 0.091 0.175 0.067 0.012 0.002 0.000 0.000 
XD 0.068 0.081 0.115 0.093 0.064 0.057 0.082 0.085 0.080 0.055 0.114 0.122 0.015 0.068 0.084 0.084 0.045 
JY 0.059 0.078 0.108 0.087 0.064 0.063 0.085 0.075 0.075 0.054 0.113 0.113 0.013 0.068 0.075 0.079 0.016 


37 茶 小 绿叶 蝉 不 同 种 群 之 间 ( 上 下 三 角 ) 和 种 群 内 部 个 体 之 间 ( 对 角 线 ) 最 近 几 代 的 基因 流 ( Nm) 1E 
Table 7 Gene flow (Nm) values over recent generations among different populations of Empoasca onukii 
( below and above diagonal) and among individuals ( on diagonal) in each population 
ЛЯ 迁 出 种 群 Source populations 


Recipient 


populations FX MB SL SG HL XH HM WZ HZ XG JD FQ JB YY IT MS XD Jy 








FX 0.790 0.007 0.007 0.007 0.007 0.126 0.007 0.007 0.007 0.127 0.007 0.007 0.009 0.007 0.007 0.007 0.008 0.011 
MB 0.011 0.687 0.007 0.007 0.019 0.007 0.007 0.008 0.008 0.007 0.008 0.008 0.007 0.007 0.007 0.007 0.007 0.007 
SL 0.007 0.008 0.728 0.014 0.007 0.007 0.051 0.007 0.007 0.007 0.007 0.007 0.007 0.028 0.009 0.088 0.007 0.007 
SG 0.007 0.008 0.059 0.680 0.007 0.007 0.014 0.007 0.007 0.007 0.007 0.007 0.007 0.027 0.010 0.087 0.007 0.007 
HL 0.015 0.180 0.007 0.007 0.854 0.007 0.007 0.007 0.007 0.007 0.007 0.007 0.008 0.011 0.007 0.008 0.007 0.007 
XH 0.048 0.009 0.007 0.007 0.009 0.715 0.007 0.007 0.007 0.055 0.007 0.007 0.008 0.007 0.007 0.007 0.007 0.008 
HM 0.007 0.008 0.054 0.057 0.007 0.007 0.725 0.007 0.007 0.007 0.007 0.007 0.007 0.027 0.013 0.089 0.007 0.007 
WZ 0.009 0.020 0.009 0.009 0.175 0.009 0.009 0.675 0.009 0.007 0.009 0.015 0.007 0.007 0.007 0.007 0.007 0.007 
HZ 0.088 0.011 0.007 0.007 0.067 0.035 0.007 0.007 0.677 0.007 0.011 0.014 0.007 0.007 0.007 0.007 0.007 0.007 
XG 0.054 0.009 0.007 0.007 0.010 0.049 0.007 0.009 0.041 0.723 0.007 0.007 0.007 0.007 0.007 0.007 0.007 0.007 
JD 0.097 0.010 0.007 0.007 0.052 0.038 0.007 0.007 0.007 0.044 0.674 0.014 0.007 0.007 0.007 0.007 0.007 0.007 
FQ 0.014 0.019 0.015 0.015 0.082 0.015 0.014 0.009 0.007 0.015 0.007 0.681 0.007 0.007 0.007 0.007 0.007 0.007 
JB 0.007 0.008 0.007 0.007 0.007 0.007 0.007 0.009 0.007 0.007 0.007 0.014 0.742 0.008 0.007 0.011 0.077 0. 108 
YY 0.015 0.008 0.014 0.013 0.008 0.009 0.024 0.009 0.007 0.010 0.007 0.014 0.007 0.699 0.035 0.033 0.007 0.007 
LT 0.007 0.008 0.009 0.013 0.008 0.007 0.022 0.009 0.007 0.007 0.007 0.014 0.007 0.029 0.701 0.030 0.007 0.007 
MS 0.011 0.008 0.012 0.013 0.011 0.008 0.021 0.009 0.007 0.009 0.008 0.014 0.008 0.123 0.139 0.799 0.007 0.008 
XD 0.008 0.008 0.007 0.007 0.007 0.007 0.007 0.009 0.007 0.007 0.007 0.015 0.040 0.007 0.007 0.007 0.709 0.112 
JY 0.007 0.008 0.007 0.007 0.007 0.007 0.007 0.009 0.007 0.007 0.007 0.015 0.077 0.007 0.007 0.008 0.042 0.779 
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R8 景观 因子 与 茶 小 绿叶 蝉 种 群 遗传 分 化 的 Simple Mantel 和 Partial Mantel 相关 性 分 析 结 果 
Table 8 Results of Simple Mantel and Partial Mantel tests among landscape factors 
and genetic differentiation of Empoasca onukii populations 
变量 分 析 类 型 Provesti 氏 遗 传 距离 Nei 氏 遗 传 距离 И 
vi uy ides Provesti' s genetic distance Nei' s genetic distance Pairwise F gp 
相关 系数 РАН P-value 相关 系数 РАН P-value 相关 系数 РАН P-value 

geo Simple Mantel 0. 2938 0.0088 ** 0.2015 0. 0455 * 0. 2052 0. 0301 * 
LCOVI Partial Mantel 0.0019 0.9893 —0.0284 0. 8401 一 0.0297 0. 8037 
LCOV2 Partial Mantel 0.0066 0.9506 —0. 0249 0. 8703 – 0. 0280 0.8104 
LCOV3 Partial Mantel 0.0003 0.9974 —0. 0367 0. 8537 -0.0395 0.7893 
elevation Partial Mantel —0. 1621 0.2139 —0. 1234 0.3128 —0. 1022 0.3562 
500m_Tea plantation Partial Mantel 0. 2568 0. 0126 * 0. 1944 0. 0404 * 0. 2700 0. 0096 ** 
500m_Cropland Partial Mantel 0.2579 0.0102 * 0.1953 0.0426 * 0. 2708 0. 0089 ** 
500m_Grassland Partial Mantel 0.2576 0.0101 * 0.1952 0. 0403 * 0.2709 0.0111* 
500m, Woodland Partial Mantel 0.2567 0. 0110 * 0. 1942 0. 0422 * 0. 2696 0. 0093 ** 
1000m_Grassland Partial Mantel 0.2576 0.0122 * 0.1953 0.0396 * 0.2709 0. 0087 ** 
1500m_Woodland Partial Mantel 0.2571 0.0084 ** 0.1946 0. 0376 * 0. 2701 0.0100 * 
2500m_Tea plantation Partial Mantel 0.2576 0.0097 ** 0.1952 0. 0385 * 0. 2708 0. 0094 ** 
2500m_noveg Partial Mantel 0.2576 0.0095 ** 0.1953 0. 0392 * 0.2709 0.0102 * 
2500m_Cropland Partial Mantel 0.2576 0.0117 * 0.1952 0. 0367 * 0.2709 0.0118 * 
2500m_Grassland Partial Mantel 0.2576 0.0102 * 0. 1953 0. 0363 * 0.2709 0.0093 ** 
2500m, Woodland Partial Mantel 0.2578 0.0121 * 0.1954 0.0401 * 0.2709 0. 0099 ** 


























*P«0.05; *P«0.01; ** P«0.001. 变量 所 对 应 的 具体 信息 见 表 3。Details for variables are listed in Table 3. 


距离 矩阵 多 元 回归 分 析 表 明 , 地 理 距 离 (gseo) 和 ”遗传 距离 Мег 氏 遗 传 距离 和 遗传 分 化 指数 的 影响 


采样 点 周转 














半径 1000 m 范围 内 草地 面积 占 比 ” 都 是 显著 的 ;这 两 个 因子 对 Provesti RHEA Nei 














(1000m_Grassland ) 两 个 景观 环境 因子 无 论 是 组 合 。 氏 遗 传 距离 和 遗传 分 化 指数 的 解释 比例 有 差异 , 决 
在 一 起 还 是 各 自 单独 对 茶 小 绿叶 蝉 种群 Provesti [E 定 系数 R 分 别 为 14.70% , 7.72% 和 12.28%( 表 9)。 











RI 距离 矩阵 多 元 回归 分 析 相 关 的 景观 因子 及 其 对 茶 小 绿叶 蝉 种 群 遗传 变异 的 解释 比例 
Table 9 Related landscape factors and their explanatory proportion to the genetic variation of 


Empoasca onukii populations by multiple regression on distance matrix analysis 








divers 回归 系数 
响应 变量 模型 解释 变量 决定 系数 R ў Е 59 РАН 
egression 
Response variables Model Explanatory variables R? P-value in F-test 
coefficient B 
HEIER 0. 1470 0.0018 ** 
Provesti 氏 遗 传 距离 
GD ~ geo + 1000m_Grassland geo 2.28E-06 0. 0063 ** 
Provesti' s genetic distance 
1000m. Grassland 5.01E-03 0.0114 * 
到 遗传 距离 0.0772 0. 0262 * 
Nei 氏 遗 传 距离 
GD ~ geo + 1000m_Grassland geo 3. 19E-06 0. 0380 * 
Nei’ s genetic distance 
1000m_Grassland 7.92Е-03 0. 0417 * 
x 0.1228 0.0069 ** 
遗传 分 化 指数 
"D GD ~ geo + 1000m  Grassland geo 8. 00E-07 0.0267 * 
Pairwise F'sr 
1000m. Grassland 2. 68E-03 0. 0083 ** 


GD; Provesti 氏 遗 传 距离 或 Nei 氏 遗 传 距离 或 遗传 分 化 指数 Provesti’ s genetic distance or Nei’ s genetic distance or pairwise Fop. 解释 变量 信息 见 
8 ST 





Æ 3, Information for explanatory variables was shown in Table3. * Р «0.05; "P «0.01; ** P «0.001. 


一 般 线性 模型 分 析 显 示 ,采样 点 周围 半径 1 000 指数 线性 关系 显著 (图 3) ,种 群 遗 传 距离 随 着 地 理 
m 范围 内 草地 面积 占 比 和 地 理 距离 与 茶 小 绿叶 蝉 种 ”距离 的 增加 而 显著 变 大 (7 = 0. 2015 ~ 0. 2938, P = 















































PE Provesti РЕНН Ej Nei 氏 遗 传 距离 和 遗传 分 化 0.0002 ~0.0120)( 图 3: A, B, C), 也 随 着 采样 点 
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与 地 理 距 离 和 采样 点 周围 半径 1 000 m 范围 内 草地 面积 占 比 的 线性 关系 


Fig. 3 








Linear correlation between Provesti' s genetic distance (A, D), Nei' s genetic distance (B, E) and pairwise For(C, F) 


of Empoasca onukii populations and geographical distance and the area percentage of grassland in the circular sector 


of 1 000 m radius around sampling site based on generalized linear mixed model 





周围 半径 1 000 m 范围 内 草地 面积 占 比 的 变 大 而 显 
著 变 大 (r=0. 1866 ~0. 2604, P 20.0015 ~0. 0226) 
(ЖЗ: D, ,下 )。 同 时 ,采样 点 周围 半径 1 000 m 
范围 内 草地 面积 占 比 与 茶 小 绿叶 蝉 种 群 He 和 Dj 之 
间 的 线性 关系 也 是 显著 的 (图 4: А, C), 茶 小 绿叶 
蝉 种 群 遗 传 多 样 性 水 平 随 着 采样 点 周围 半径 1 000 
m 范围 内 草地 面积 占 比 的 变 大 而 变 大 (7 = 0. 4829 ~ 












































0. 5222, P 20.0267 ~0.0426)。 
3 讨论 
3.1 茶 小 绿叶 蝉 的 自主 扩散 能 


本 研究 基于 23 对 微 卫星 位 点 的 基因 型 数据 对 
福建 省 安溪 县 18 个 茶 小 绿叶 蝉 种 群 进行 了 分 析 , 虽 




















10} 李 金 玉 等 : 茶园 周边 景观 格局 对 茶 小 绿叶 蝉 种 群 遗传 结构 的 影响 
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4 一 般 线性 模型 分 析 茶 小 绿叶 蝉 种 群 He (А), Ho (В), Dw(C) 和 4An(D) 与 采样 点 














周围 半径 1 000 m 内 草地 面积 占 比 的 线性 关系 


Fig. 4 Linear correlation between He (A), Ho (B), D,,(C) and A,(D) of Empoasca onukii populations and the area percentage 

















of grassland in the circular sector of 1 000 m radius around sampling site based on generalized linear mixed model 


Не: 期 望 杂 合 度 Expected heterozygosity; Но: 观察 杂 合 度 Observed heterozygosity; уу: 基因 多 样 性 Gene diversity; Ap: 等 位 基因 丰富 度 Allelic 


richness. 


然 是 小 尺度 区 域 的 研究 ,基于 模型 假设 和 不 基于 模 
型 假设 的 分 析 结 果 显 示 研 究 区 域内 茶 小 绿叶 蝉 种 群 
遗传 分 化 结构 明显 ,遗传 距离 和 分 化 指数 的 分 析 显 
示 大 部 分 种 群 间 的 遗传 分 化 达到 了 中 等 水 平 ,多 元 
回归 模型 和 一 般 线性 模型 以 及 Mantel 分 析 显 示 茶 
小 绿叶 蝉 种群 遗传 分 化 与 地 理 距 离 显 著 正 相关 ,种 
群 表 现 地 理 距 离 隔 离 格 局 ,这 些 信息 表 明 茶 小 绿叶 
蝉 种 群 有 明显 的 遗传 分 化 格局 。 近 期 有 其 他 研究 者 
报道 了 基于 微 卫 星 数据 分 析 茶 小 绿叶 蝉 在 中 国 四 大 
茶 区 明显 的 种 群 遗 传 结构 和 地 理 距 离 隔离 格局 
( Zhang et al., 2019) , Lo 等 (2019 ) 基于 微 卫星 标记 
发 现 澳 大 利 亚 16 МЛН Thaumastocoris peregrinus 
种 群 表现 中 等 分 化 和 地 理 距 离 隔 离 格局 并 据 此 推测 
枸 椿 的 自主 扩散 能 力 有 限 ;Cao 2 (2019 ) 对 15 个 瓜 
fij CL Thrips palmi 种 群 的 微 卫 星 基因 型 分 析 发 现 种 
群 表 现 出 显著 的 遗传 分 化 和 地 理 距离 隔离 格局 , 基 









































































































































于 BAYESASS 估算 的 种 群 间 的 现时 基因 流明 显 受 
限 ,推测 种 群 间 可 能 没有 持续 的 扩散 事件 。 凡 是 具 
有 远 距离 和 大 范围 扩散 的 迁 飞 性 昆虫 种 群 ,往往 不 
具有 明显 的 遗传 结构 ,异地 种 群 间 通常 可 以 检测 到 
明显 的 基因 流 。 例 如 , 迁 飞 性 害虫 小 菜 蛾 Plutella 
xylostella 种 群 的 遗传 分 化 不 明显 ,具有 明显 的 基因 
流 格局 ( Wei et al., 2013, 2017) ;而 采 自 相同 区 域 的 
HMPA A AWE Cotesia vestalis 的 种 群 间 却 表现 出 很 
强 的 遗传 分 化 水 平 (Wei et al., 2017) ,因此 认为 菜 
峨 盘 绒 草 峰 可 能 不 会 跟随 其 寄主 远 距 离 迁 飞 。 
Zhang 等 (2019 ) 虽然 未 对 基因 流 进行 直接 的 佑 算 与 
分 析 , 但 认为 茶 小 绿叶 蝉 不 同 地 理 种 群 之 间 存 在 大 
量 基因 流 。 本 研究 基于 BAYESASS Xp np tg 
异地 种 群 间 的 迁移 率 ( 基 因 流 ) 进行 了 估算 ， 
BAYESASS 估算 的 迁移 率 (m) 通 常 被 认为 是 种 群 5 
代 以 内 (<5 代 ) 的 情况 (Wilson and Rannala, 



































































































































1256 昆虫 学 报 Acta Entomologica Sinica 63 35 





2003) ,鉴于 茶 小 绿叶 蝉 种 群 在 中 国 茶 区 每 年 发 生 
9 - 16 代 (Fu et al., 2014) ,本 人 研究 估算 的 基因 流 可 
以 反映 茶 小 绿叶 蝉 种 群 “当代 ”( 不 到 一 年 ) 的 迁移 
扩散 事件 ,也 即 帮助 我 们 推测 茶 小 绿叶 蝉 种 群 的 扩 














定 树龄 的 茶园 采集 样本 以 避免 近期 茶 苗 移植 对 分 析 
的 影响 ,同时 选用 微 卫星 分 子 标记 和 BAYESASS 分 
析 种 群 遗 传 信息 。 一 般 认 为 高 度 可 变 的 位 点 ,如 微 
卫星 分 子 标记 适合 于 研究 和 反映 生物 种 群 短期 的 小 














散 过 程 和 迁移 能 力 。BAYESASS 分 析 显 示 其 研究 区 
域内 茶 小 绿叶 蝉 种 群 间 的 当代 基因 交流 受 限 ,其 现 
时 基因 流 往往 仅 限 于 种 群 内 部 个 体 之 间 , 表 明 小 尺 
度 区 域 范围 内 茶 小 绿叶 蝉 异 地 种 群 间 可 能 缺乏 自主 
的 扩散 和 交流 过 程 。 根 据 以 上 信息 ,我们 推测 茶 小 
绿叶 蝉 很 可 能 不 具有 远 距离 迁 飞 的 习性 ,自主 扩散 
的 能 力 有 限 。 周 宁 宁 (2014) 和 陈 世 春 等 (2015 ) 分 
别 基于 线粒体 COI 和 16S rRNA 基因 序列 对 假 眼 小 
绿叶 暗 Empoasca vitis 地 理 种 群 的 分 析 也 都 认为 茶 
小 绿叶 蝉 可 能 是 一 类 扩散 能 力 有 限 的 害虫。 
尽管 茶 小 绿叶 暗 的 自主 扩散 能 力 有 限 , 人 类 活 
动 对 茶 小 绿叶 蝉 种 群 的 影响 是 不 可 忽视 的 ,其 中 最 
主要 的 是 茶 苗 在 异地 间 的 调运 和 移植 过 程 可 能 可 以 
为 茶 小 绿叶 蝉 种 群 实现 跨 地 区 扩散 和 远 距 离 传 播 创 
造 很 好 的 条 件 和 机 会 。 在 过 去 几 十 年 中 , 随 着 无 性 
繁殖 技术 在 茶树 上 的 运用 和 发 展 ,由 于 其 相对 于 茶 
籽 播 种 来 说 更 加 快速 高 效 , 越 来 越 多 新 茶园 的 建立 
选择 从 异地 调运 茶 苗 ,这 为 茶 小 绿叶 蝉 种 群 实现 跨 
地 区 扩散 和 远 距 离 传 播 创造 了 很 好 的 条 件 和 机 会 。 
茶 小 绿叶 蝉 种 群 具 有 跟随 寄主 茶 苗 实现 跨 地 区 扩散 
和 远 距 离 传 播 的 能 力 , 茶 小 绿叶 蝉 通常 将 卵 产 入 茶 
术 嫩 梢 或 叶 梢 组 织 中 ,其 成 幢 和 知 虫 也 很 小 (Shi et 
al., 2015) ,又 具有 和 久 坐 不 动 和 大 部 分 时 间 均 栖息 在 
叶 背 或 嫩 梢 上 吸食 汁液 为 害 的 习性 ( 边 夭 等 ， 
2014) 。 在 种 植 者 没有 更 多 农事 操作 和 育苗 知识 的 
情况 下 , 茶 苗 很 可 能 受到 茶 小 绿叶 蝉 侵 染 ,并 使 其 跟 
随 茶 苗 异 地 调运 而 实现 远 距 离 传 播 和 异地 扩散 。 
此 ,新 建 茶园 仍 要 关注 茶 小 绿叶 蝉 侵 染 和 危害 的 问 
题 。 随 着 未 来 新 茶园 的 发 展 和 旧 茶 园 的 蔡 换 , 茶 苗 
和 茶树 在 异地 间 的 调运 应 进行 严格 的 检查 以 防止 茶 
小 绿叶 蝉 种 群 在 更 大 区 域 范围 内 的 扩散 和 传播 。 同 
时 , 相 比 于 择 捕 等 无 性 繁殖 技术 ,组 培 和 茶 籽 播种 可 
以 繁育 出 不 带 虫 卵 的 茶 苗 ,结合 将 这 些 “ 无 虫 苗木 
在 温室 中 培育 的 措施 可 以 帮助 减少 潜在 的 茶 小 绿叶 
蝉 侵 染 的 风险 ,并 避免 更 大 范围 的 传播 和 扩散 。 
虽然 人 类 活动 影响 不 可 忽视 也 不 可 控 , 但 茶树 
是 多 年 生 作物 ,一 经 种 植 可 以 多 年 保持 ,因此 在 已 经 
建立 的 成 熟 茶园 中 ,生态 系统 相对 于 其 他 作物 系统 
更 稳定 ,人 类 活动 对 茶 小 绿叶 蝉 的 助 迁 作用 和 机 会 
会 大 大 减少 。 因 此 研究 中 我 们 选择 其 中 茶树 具有 一 











































































































尺度 区 域 的 遗传 结构 和 基因 流 ( Wang, 2010, 
2011)。 而 基于 频率 分 配 检测 的 方法 BAYESASS 被 
广泛 用 于 估算 当前 或 以 前 几 代 的 迁移 率 来 推断 最 近 
的 基因 流动 。 由 此 研究 中 分 析 的 当代 种 群 遗传 分 化 
结构 和 基因 流 空 间 格局 主要 反映 茶 小 绿叶 蝉 种 群 的 
自主 扩散 或 迁 飞 能 力 信息 ,而 可 能 与 茶 苗 异地 调运 
过 程 关系 较 少 。 

3.2 茶园 周边 景观 格局 对 茶 小 绿叶 蝉 种 群 遗传 结 
构 的 影响 

本 研究 中 ,研究 区 域内 茶 小 绿叶 蝉 种 群 大 体 上 
分 为 东部 西北 部 和 南部 种 群 的 格局 ,结合 茶园 周边 
景观 格局 分 析 结果 和 采样 观察 记录 以 及 相关 文献 记 
录 ( 税 伟 等 , 2017) ,位 于 南部 的 JB, LT 和 YY 等 种 
群 所 在 乡镇 ,主要 采用 单一 化 和 和 集约 化 的 种 植 模式 ， 
在 茶园 周围 少见 自然 和 半 自 然 生 境 或 其 他 作物 田 ; 
位 于 西北 部 的 FQ 和 XG 等 种 群 所 在 乡镇 ,茶园 大 都 
分 布 于 高 海拔 的 山坡 上 ,周围 经 常会 有 大 量 森 林 等 
自然 和 半 自 然 生境 分 布 。 而 位 于 东部 的 FX, SG, 
HL, SL 和 MB 等 种 群 所 在 的 乡镇 ,种 植 模式 则 会 更 
多 样 化 ,茶园 周围 经 常会 有 更 多 居民 点 和 其 他 作物 
田 的 分 布 。 根 据 这 些 信息 ,我 们 推测 周边 景观 格局 
不 同 的 茶园 中 的 茶 小 绿叶 蝉 种 群 遗 传 结构 和 基因 交 
流 模式 可 能 存在 差异 。 

丝 助 于 景观 遗传 学 的 原理 和 方法 ,使 用 Mantel 
相关 性 分 析 、 距 离 矩 阵 多 元 回归 模型 MRM 和 一 般 
线性 模型 GCLMM 定量 化 解析 了 地 理 距 离 海拔、 样 
点 周围 不 同 尺度 范围 的 景观 组 成 因子 以 及 景观 结构 
与 茶 小 绿叶 蝉 种 群 遗传 变异 的 相关 性 。 结 果 显 示 地 
理 距离 (geo) 和 采样 点 周围 半径 1 000 m 范围 内 草 
地 面积 占 比 (1000m_Grassland ) 是 与 研究 区 域内 茶 
小 绿叶 蝉 种 群 遗 传 分 化 显著 相关 的 因子 , 茶 小 绿叶 
蝉 种 群 遗传 分 化 水 平 随 着 采样 点 间 地 理 距 离 的 变 大 
而 变 大 ,也 随 着 采样 点 间 周 围 草地 面积 差异 值 的 变 
大 而 显著 变 大 ,与 周围 草地 面积 较 小 的 茶园 相 比 ,有 
较 多 草地 分 布 的 茶园 中 的 茶 小 绿叶 蝉 种 群 表现 出 较 
高 的 遗传 多 样 性 水 平 。 类 似 地 ,Dong 等 (2018 ) 的 研 
SUA PU HS Aphis gossypii 的 遗传 多 样 性 水 平 与 
1 000 m 和 2 000 m 范围 内 的 草地 面积 比例 呈 显 著 
正 相 关 。 已 有 的 大 量 景观 生态 学 研究 表明 一 些 生 境 
(special habitat) 或 土地 利用 类 型 (land cover) 会 影 
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响 昆 虫 的 定 殖 过 程 (Dong et al., 2015; Angelella et 
al., 2016) 。 因 此 ,景观 组 成 会 影响 昆虫 种 群 的 建立 
和 扩张 过 程 ,进而 也 可 能 影响 其 种 群 的 遗传 结构 和 
分 化 (Dong et al., 2018) ,特别 是 对 于 这 类 自主 扩散 
能 力 有 限 的 昆虫 来 说 ,景观 组 成 的 影响 作用 可 能 会 
更 明显 。 




















那 杂 草地 生境 对 茶 小 绿叶 蝉 种 群 的 促进 作用 和 在 保 
育 天 政 资 源 以 更 好 发 挥 防 控 效 能 作用 间 的 比较 和 权 
衡 也 是 一 个 值得 探讨 的 复杂 问题 。 

此 外 ,景观 尺度 可 能 也 对 茶 小 绿叶 蝉 种 群 的 扩 
散 和 定 殖 有 一 定 的 影响 ,本 研究 在 5 个 空间 尺度 
(景观 半径 0.5 ~2.5 km 不 等 ) 上 研究 了 景观 组 成 














本 研究 中 的 土地 利用 类 型 空间 分 布 格局 可 以 简 
单反 映 茶 小 绿叶 蝉 的 栖 境 和 寄主 植物 以 及 不 适宜 生 
境 的 分 布 。 已 有 的 研究 还 不 能 排除 茶 小 绿叶 蝉 会 在 


因子 对 茶 小 绿叶 蝉 种 群 遗传 多 样 性 和 变异 的 影响 ， 
结果 显示 茶园 周围 1 000 m 半径 内 的 草地 面积 解释 
了 茶 小 绿叶 蝉 种 群 遗传 多 样 性 和 变异 的 最 大 比例 。 





























杂 草 上 越冬 .过渡 和 繁殖 的 可 能 性 ( 陈 李 林 等 ， 
2019; Chen et al., 2019)。 林 草地 对 茶 小 绿叶 暗 种 
群 的 遗传 多 样 性 具有 积极 作用 ,可 能 一 方面 茶 小 绿 
叶 暗 种 群 在 茶园 受到 农事 操作 ( 如 茶 校 修剪 和 农药 
喷 酒 等) 干扰 时 移动 到 杂 草 地 过 渡 , 干 扰 解 除 后 再 
从 杂 草 地 移 回 茶园 栖息 地 ,景观 组 成 可 能 通过 这 样 
的 方式 影响 其 种 群 的 建立 、 定 殖 过 程 ,进而 影响 其 遗 
传 结构 和 分 化 ; 另 一 方面 , 杂 草 地 生境 本 身 也 可 能 提 
供 丰 富 的 叶 暗 资源 ,寄主 和 替代 寄主 分 布 通过 选择 
和 繁殖 有 影响 茶 小 绿叶 蝉 种 群 的 遗传 结构 和 分 化 。 相 
比 于 杂 草 地 ,森林 生境 包含 更 多 的 多 年 生 乔 木 树种 ， 
这 些 高 大 的 树种 可 能 大 多 数 都 不 能 成 为 茶 小 绿叶 暗 
合适 的 寄主 植物 ,而 且 随 着 森林 面积 的 增加 ,适合 于 
茶 小 绿叶 蝉 种 群 栖息 的 茶园 面积 减少 ,同时 也 不 能 
排除 森林 生境 中 的 高 大 树种 会 阻 断 茶 小 绿叶 蝉 种群 
间 的 交流 .扩散 和 连续 性 的 可 能 性 。 因 此 ,在 周围 有 
更 多 条 林 生 境 分 布 的 样 点 中 的 种 群 如 FQ 和 XG SF, 
叶 暗 种 群 表现 较 低 遗传 多 样 性 水 平 ,一 般 线 性 模型 
的 分 析 也 显示 叶 暗 种群 遗传 多 样 性 水 平 与 茶园 周转 
的 森林 生境 面积 呈 负 相关 (数据 未 展示 )。 赵 紫 华 
等 (2010) 报 道 与 麦田 相 邻 的 自然 和 半 自 然 生境 枉 
境 可 有 效 地 延缓 麦 蚜 种 群 迁 入 作物 田 的 时 间 并 减少 
扩散 传播 事件 。 从 本 研究 的 结果 来 看 ,不 同 的 自然 
生境 类 型 对 害虫 种 群 的 影响 作用 有 差异 。 就 茶 小 绿 
叶 蝉 而 言 ,茶园 周围 的 杂 草 地 生境 对 其 种 群 建立 .过 
渡 和 扩张 可 能 是 有 利 的 ,因此 清除 杂 草 可 能 是 茶 小 
绿叶 蝉 可 持续 治理 策略 很 重要 的 一 项 内 容 。 当 然 ， 
非 作物 生境 通常 也 被 认为 可 能 栖息 着 大 量 可 作为 生 
物 防 治 资源 的 天 敌 物 种 (Tscharntke et al., 2012; 
Gurr et al., 2016) , 茶 小 绿叶 蝉 有 两 种 主要 的 卵 寄生 
Ws IH} du — fas | Е Stethynium empoascae 和 微小 裂 
骨 继 小 蜂 Schizophragma parvulas 在 田间 的 寄生 率 超 
iX 30% ,被 认为 是 茶 小 绿叶 蝉 可 持续 控制 的 重要 天 
T(E AMI, 2008; 毛 迎 新 等 , 2008) ,那么 
Aem ENDE X PER AE NERO aT 如 果 有 益 ， 
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尽管 当前 我 们 对 茶 小 绿叶 蝉 种 群 的 移动 特性 还 了 解 
得 不 多 ,但 茶园 周围 1 000 m 半径 这 个 尺度 范围 的 
景观 可 能 包含 了 茶 小 绿叶 蝉 种 群 的 生态 邻 域 
(ecological neighborhood) ( Dong et al., 2018) , 其 中 
可 能 包含 或 者 满足 了 其 可 以 或 需要 利用 的 生境 
RE. 

综合 以 上 信息 来 看 ,关于 小 尺度 区 域 范围 内 茶 
小 绿叶 蝉 景观 遗传 学 的 研究 可 能 有 助 于 增进 对 茶 小 
绿叶 蝉 定 殖 及 其 相应 空间 尺度 和 相关 景观 格局 的 理 
解 。 但 是 ,生物 种 群 遗 传 变 异 是 一 个 复杂 的 过 程 , 微 
卫星 分 子 标记 ( SSR ) 的 统计 功效 也 还 比较 有 限 , 研 
究 工 作 还 比较 初步 ,16 个 空间 因子 之 外 ,可 能 还 有 
其 他 相关 的 影响 因子 没有 被 检测 或 识别 到 。 未 来 的 
工作 可 以 考虑 对 茶 小 绿叶 蝉 种 群 及 其 相关 天 敌 昆 虫 
种 群 在 茶园 景观 中 的 生态 学 过 程 进行 更 多 的 观察 和 
研究 ,如 检查 可 能 充当 茶 小 绿叶 蝉 替 代 寄 主 的 植物 
物种 测量 从 茶园 栖 境 迁 出 和 从 过 渡 栖 境 迁 回 茶园 
的 最 近 ( 远 ) 距 离 等 等 ,以 进一步 确认 茶园 周边 景观 
格局 对 茶 小 绿叶 蝉 种 群 建立 、 定 殖 、 扩 张 的 影响 。 
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